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Abstract 

In this paper we represent a new form of condition for the consistency of the 
matrix equation AXB = C. If the matrix equation AXB = C is consistent, 
we determine a form of general solution which contains both reproductive and 
non-reproductive solutions. Also, we consider applications of the concept of 
rcproductivity for obtaining general solutions of some matrix systems. 
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1. Introduction 

The concept of the reproductive equations was introduced by S.B. Presic [1] 
in 1968. In this part of the paper we give the definition of reproductive equations 
and the most important statements related to the reproductive equations. Using 
the concept of reproductivity in the next section we obtain the general solutions 
of some matrix systems. 

Let S be a given non-empty set and J be a given unary relation of S. Then 
an equation J(x) is consistent if there is at least one element Xq £ S, so-called 
the solution, such that J(xq) is true. A formula x — </>(t), where <f> : S — ¥ S is a 
given function, represents the general solution of the equation J(x) if and only if 

(Vt)J(<j>{t)) A (Vaj)(J(x) => {3t)x = 4>{t)). 

In this part of the paper we give the definition of reproductive equations and 
the fundamental statements related to the reproductive equations. 



2010 Mathematics Subject Classification: 15A24 
Email addresses 

Corresponding author : Biljana Radicic <radicic.biljana@yahoo.com> ( "part time job), 
Branko Malesevic <malcsevic@etf.rs> 



Preprint submitted to arXiv.org 



May 13, 2012 



Definition 1.1. The reproductive equations are the equations of the following 

form: , s 

x = <p(x), 

where x is a unknown, S is a given set and tp : S — > S is a given function 
which satisfies the following condition: 

tpOtp = Lp. (1) 

The condition ((T|) is called the condition of reproductivity. The fundamental 
properties of the reproductive equations are given by the following two state- 
ments S.B. Presic Q (see also [Bj], [H and 0). 

Theorem 1.1. For any consistent equation J(x) there is an equation of the 
form x — tp{x), which is equivalent to J{x) being in the same time reproductive 
as well. ♦ 

Theorem 1.2. If a certain equation J{x) is equivalent to the reproductive one 
x = <p(x), the general solution is given by the formula x = (f(y), for any value 
yeS.* 

Let us remark that a formula x — <p(t) , where (j) : S — > S is a given function, 
represents the reproductive general solution [29j j of the equation J{x) if and only 

(W)J(#*)) A (W)(J(t) =*•* = #*)). 

S.B. Presic was the first one who considered implementations of reproductiv- 
ity on some matrix equations 0] (see also @], @)- The concept of reproductivity 
allows us to analyse various forms of the solution. General applications of the 
concept of reproductivity were also considered by J.D. Keckic in Q, [H and J.D. 
Keckic and S.B. Presic in [lq ]. 



2. The matrix equation AXB — C 



Let m,n GN and C is the field of complex numbers. The set of all matrices 
of order m x n over C is denoted by C mx ™. We use denotement C™ x ™ for the 



set of all m x n complex matrices whose rank is a. Let A = [c. 



\ By 



A T , Ai^ and A^j we denote the transpose of A, the i-th row of A and the j-th 
column of A, respectively. 
Therefore, 

and 



A 
A 



i.i 



= (ai,i,a ii2 , Oi jn ), i = 1, .. 

a 2,j, «raj) T , j = 1) 



A solution of the matrix equation AX A = A is called {l}-inverse of A and it is 
denoted by A^\ 
The matrix equation 

AXB = C (2) 

was considered by many authors (see [1 3|. [14| . 
many papers (see [TO], 0, 53, HTHJ-Sc 



25]. 33]). There are 



341) where the matrix 



equation @ is studied as a part of different matrix systems. 
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This part of the paper is organized as follows: In 2.1. we represent the new 
form of condition for the consistency of the matrix equation ([2]) . An extension 
of Penrose's theorem related to the general solution of the matrix equation © 
is given in 2.2. Namely, we represent the formula of general solution of the 
matrix equation © if any particular solution Xq is known. In 2.3. we give 
a new form of particular solution Xq, using the result in 2.1., such that the 
formula of general solution of the matrix equation ([2]), which is given in 2.2., is 
reproductive. In 2.4- we give two applications on some matrix systems which 
are in relation to the matrix equation ©. 



2.1. Let A e C™ x '\ B e €-1 q and C £ C mxq . The matrix A E C r " x ™ has a 
linearly independent rows. If the linearly independent rows of the matrix A 
are not at the first a positions, by multiplying the matrix A by a permutation 
matrix Ta on the left, we can permute the rows of the matrix A. 

Therefore, the matrix 

A = T A A (3) 
has linearly independent rows at the first a positions. Analogue, the matrix 

B = BT B (4) 
has linearly independent columns at the first b positions. 

The considerations which follow are valid for any choice of matrices Ta and Tb 
such that A has linearly independent rows at the first a positions and B has 
linearly independent columns at the first b positions. 

Let 

C = T A CT B . (5) 

Next, let for the matrices A and B regular matrices Qi,Pi,Q 2 and Pi be de- 
termined such that the following equalities are true: 



Q 1 AP 1 = E a 





" 









and 



Q 2 BP 2 =E b = 



' h 


" 









(G) 



i.e. 

A = Q^EaP^ 1 and 
Then, from ©, Q and © we get that 

A = TaQi 1 EaP^ 1 and 



P — Ql 1 PbP 2 • 



p — Q 2 l PbP 2 1 Tb 



(7) 



i.e. 



A=(Q 1 T A 1 )- 1 E a P^ 1 and 
If we introduce the following designations: 
Ql = QiTJ 1 and 



P = Q 2 1 E b (T B 1 P 2 )-\ 



P2 = Tg 1 P 2 



(8) 



3 



we get that 

A = Ql~ 1 E a P^ 1 and B = Q^ 1 E b P 2 ~\ (9) 
Considering that the following equalities are true for A^> and B^> : 



A<U = P l 



la 


Xx ' 


x 2 


x 3 



Qx and = P 2 



lb 




Y 2 


Y 3 . 



(10) 



where X\, X 2l X 3 , Yx, Y 2 and Y3 are arbitrary matrices corresponding dimen- 
sions (see C. Rohde [3j), then: 



la 


Xx ' 


x 2 


x 3 



i.e. 



and 



1« = Px 



AA^=Qx EaP^Px 



B^B = P 2 



Qx and B {1) = % 



Qi = Qi 



h 


Yx ' 


Y 2 


Y 3 . 



Q2 



la 


Xx ' 


x 2 


x z 



h 


Yx " 


Y 2 


Y 3 . 



hQ^EbK = P~ 2 



' la 


Xx 








" h 


" 


I ^ 






Pa 



(11) 



(12) 



As we mentioned, the matrix A has linearly independent rows at the first a 
positions and the matrix B has linearly independent columns at the first b 
positions. 

Let 



Ai^ = ^ a i,i A i 



and 



b 



Bij ■= j = b+l,...,q. 



(13) 
(14) 



k=l 



for some scalars a^; and ftkj- Then, Qx and P 2 have the following form: 





Qi = 


' la 







and P 2 = 


' h 


L 2 


(15) 




[ Lx 











Iq-b 


where 






















— OL a +xA 


.. -a a 


+1,0 








~Pl,b+l 


- ~Pi, q ' 


Lx = 










and L 2 = 














,1 


—a 


m,a 








-flb,b+l 


- _ 
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Let us emphasize that the following statement is true. 



Lemma 2.1. Let AG< 



5eCp, CeC mxq . Suppose that A and B 



determined by (3) and (4). Then, the conditions 

AA^CB^B = C 



and 

are equivalent. 



AA^CB^B = C 



(16) 
(17) 



Proof. AA^CB^B = C ^=4> T^ l AA^T A CT B B^ BT B l = C 
AA^T A CT B B {l) B = T A CT B ^ AA^CB^B = C. ♦ 

In the following statement we give the necessary and sufficient condition so that 
the condition (fTBf is true for any choice of {l}-inverses AW and B^\ 

Theorem 2.1. Let A 6 C™ xn } B e C p b xq , C € C mX9 . Suppose that A and B 
are determined by {3p and ^ and that M3\) and |j^[ ) are satisfied. Then, the 
condition HI 6]) is true for any choice of {1} -inverses AW and 5« iff 



G-- 



Cl,l 



Ca,l 

L = l 



Cl,6 



E Pk,b+lCl,k 



Ca.b Y.Pk.b+lC a ,k 



E a mJ c !,t E E a m,lPk,b+\C-l,k 



E Pk,qCl,k 



E Pk,qC a ,k 

b k = 1 

E E Ola+l,lPk,qCl,k 



E T. a m,lPk,qCl r 



where a j are arbitrary elements of C. 



Proof. (=^) : Suppose that the condition (TT6")) is valid for any choice of {1}- 
inverses A^ 1 ) and i^ 1 ). Based on Lemma T2. II the condition (TT71) is also valid. 
Then, considering the equalities (fTT|) and (fT2j) . we get the following equality 



1 


' /a 













Q1CP2 



" 4 


' 








P2 



c. 



By multiplying the previous equality by Qi on the left and by P2 on the right 
we get 



la 


x l ' 









QICP2 



' /ft 


' 








= Q1CP2. 



(18) 
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Suppose that 



C 



Cl, 



We are going to show that C has the following form: 
Ci,i ••• Ci jfc H@k,b+lCl,k 



c= 



Ca.l 

1 = 1 



i = l 



Let 



,b+lC a ,k 

fe=l 



From (fT5|) we obtain that 
for i = l,...,a,j = 1,...,6 
for i = l,...,a,j = b+l,...,q 

for i — a + 1, ...,m, j = 1, 6 



i=lfe= 


i 




F = 


' /a 















4 











£ Pk,qCl,k 



£ fik,qC a ,k 

fe=l 

a b 

£ E Cta+l,lPk,qCl,k 



£ E OC m ,lPk,qCl,k 

(19) 



k,j&i,kl 



for i = a + 1, m, j = 6 + 1, ...,<? = Ci,j - 2J a M c 'j 



, CtilClk) 



and 

for i = 1, a, j = 1, b 



Fi,j — Ci,j 



-^h,j{ c *,k -^2 

k=l 1 = 1 

rn a 

=a+l '=1 

q b 

+ J2 yk, 3 \- C i,k-J2^k,k c i,k 

k=b+l k=1 
m a 

+ x i jHk~Y. a n c i,k 

~l=a+l 1=1 
b a 

~ Y Pk$(%k ~ 2 a ~l,l c l,k)}] 



k=l 



1 = 1 



G 



for i = 1, a,j = b + l, q Fij = 0, 

for i = a + 1, m, j = 1, b F^j = 0, 

for i = a + 1, m, j = b + 1, q Fi.j = 0. 

Finally, from JTH]) and flU]) i.e. E — F we get that 

for i = 1, a, j = 1, b c%j are arbitrary elements of C, 

for i = a + 1, m, j = 



6 



for i = l,...,a,j = 6+ !,...,<? c,j- = ^AjC,^, 



for i = a+ 1, ...,m, j = 6+1, c hJ = aj,iPk,jCi,k- 

( < S=) : Suppose that 



;=1 fc=l 



c= 



Cl,l 



Ca,l 



E a a+l,J C i,l 

i = l 



Ca,b 



tPk 



,fe+lC a ,fe 



T. a a+l,lCl.b T. Y. a a+l,lPk.b+lCl.k 



E a m/M ■■■ Y. a m,lCl.b jZT. a m,lPk,b+lCl.k ■■■ T. Y. a m,lPk,qCl.k 



E Pk,qCl,k 
T.Pk,qC a ,k 

k = l 

E E Ot a +l,lPk,qCl.k 



Then. 



and 



" 4 


*1 " 









Q1CF2 



4 


' 


^2 






Cl,l 


■■■ Ci,6 


... 


■ 


Ca,l 


■■■ Ca.6 


6 


6 





... 


... 





6 


I 6 


6 


6 . 




Cl,l ... 


Cl,b 







Co,l ■•■ 


C a ,b 







... 










... 









(20) 



(21) 
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From and ([2TJ we conclude that 



la 


*1 " 









Q1CP2 



h 


" 


Y 2 






= Q1CP2 



By multiplying the previous equality by Q\ on the left and by P 2 on the 



right we obtain the following equality: 



Qi 



1 


' la 


x 1 ' 











Q1CP2 



' lb 


' 








= c. 



(22) 



From (|22|) , considering the equalities (|TT|) and (fT2)l . we see that the condition 
(fl7|) is true. Finally, based on Lemma [2~T1 we conclude that the condition (|16|1 
is true.4 

Example 2.1. Let be given the following matrices: 



id B = 
















1 


-3 


2 


A = 


2 


1 


-1 




-1 


-4 


3 




3 


-2 


1 






1 


2 


00' 


-1 





3 


1 


4 


2 





4 


1 


5 


3 





2 


3 


5 


-1 



Then, rank(A)=2, rank(B)=2 and for 



T A 






1 








" 




' 





1 














1 










1 























1 





and Tb = 





1 























1 













1 





1 




























1 



1 


-3 


2 


2 


1 


-1 


1 


-4 


3 


00 


-2 


1 












we get that 



T A A 



Therefore, 

A 3 ^ = - A 2 
and 

B i3 = 0B U + 0B l2 
From this we get that 

"3,1 = 1, "3,2 = - 



B = BT B = ... = 



1 


2 


3 


-1 


3 


1 


4 


2 


4 


1 


5 


3 


2 





5 


-1 



A 2 



A 5 ^ = OAl 



(L4, 



) B^4 — B^i + B±2, B± 5 — B^i — B± 2 - 
1, a 4 1 = 1, a 4i2 = 1, as i = 0, a 5 . 2 = 



8 



and 



01,3 = 0, 02,3 = 0, 01,4 = 1, 02,4 = 1, 01,5 = 1, 02,5 = -1. 



Based on Theorem \2.1\ each matrix C which has the following form 



C = 



Cl,l Ci_ 2 

02,1 C 2 ,2 

Cl,l - C2,l Ci, 2 - C 2 ,2 

01,1 + 02,1 Cl,2 + C 2 ,2 





01,1 + Cl.2 

02,1 + 02,2 
Cl,l - 02,1 + Cl.2 - C 2 ,2 
Cl,l + 02,1 + Cl.2 + C 2 .2 





Cl,l — Cl.2 

02,1 — C2.2 
Cl.l — C2,l - 01,2 + C 2 ,2 
Cl.l + 02,1 — Cl,2 - C 2 ,2 




satisfies the condition |-?7| ). From that we conclude that each matrix C 



T^CTg 1 which has the following form 



C 





Cl,l Cl,2 01,1 + Cl.2 

02,1 C 2 ,2 C2,l + C 2 ,2 

01,1 — 02,1 Ci j2 - C 2 ,2 01,1 — C2,l + Cl,2 — C2,2 

Cl,l+C2,l Cl,2 + C 2 ,2 01,1+02,1+01.2+02.2 





01,1 — Cl.2 

02,1 — C 2 ,2 
01,1 - C 2 ,l - Ci j2 + C 2 ,2 
Cl.l + C 2 ,l — 01,2 - C 2 ,2 



satisfies the condition 116]) . ♦ 



In the following example we will show one useful application of Theorem l2.ll 
Namely, the result in Theorem 1 2. II can be applied to determine the consistency 
of the matrix equation ([2]). 

Example 2.2. Let A and B are the matrices as in Examvle \2.1\ and 



a) C 





1 

2 

-2 1 

2 1 




1 

-4 
5 

-3 



b) C 





1 

2 

-2 1 

2 1 




1 

4 
5 
-3 



// we compare the matrix C from a) and from b) with the general form of the 
matrix C which satisfies the condition i ll 6]) (see Examvle \2.1\ ) we see that the 
matrix C from a) satisfies the condition 116]) and the matrix C from b) does not 
satisfy the condition HI 6}) . Therefore, the matrix equation (0) is consistent for 
the matrix C from a) and it is not consistent for the matrix C from b).^ 

2.2. Recall that the matrix equation AXB = C is marked with @. In the paper 
[1] R. Penrose proved the following theorem related to the matrix equation @. 

Theorem 2.2. The matrix equation ^ is consistent iff for some choice of {1}- 
inverses and B^ of the matrices A and B the condition Mb}) is true. The 
general solution of the matrix equation ^ is given by the formula 

X = f(Y) = AWCBW +Y- A^AYBB {1 \ (23) 

where Y is an arbitrary matrix corresponding dimensions.^ 



9 



Remark 2.1. If the matrix equation fl[) is consistent, the equivalence 

AXB = C X = f(X) = X - A (1) (AXB - C)B (1) (24) 

is true. Therefore, the starting equation is equivalent to some reproductive equa- 
tion. Based on Theorem \1.2\ we can also conclude that \23\l is the general 
solution of the matrix equation (0).O 

In this paper we give an extension of Theorem 12.21 

Theorem 2.3. If Xq is any particular solution of the matrix equation (0), the 
general solution of the matrix equation (0) is given by the formula 

X = g(Y) = X + Y - A W AYBB {1 \ (25) 

where Y is an arbitrary matrix corresponding dimensions. The function g sat- 
isfies the condition of reproductivity {Ip iff Xq — A^CB^K 



Proof. It is easily to see that the solution of the matrix equation is given 
by ((25)) . On the contrary, let X is any solution of the matrix equation ([2]), then 

x = x-aWcbW + a<ucbW 

= X - A^AXBB^ + AWAXqBBW 

= X-A^A(X-X )BB^ 

= X a + (X - X ) - A {1) A(X - X )BB {1) 

= Xq + Y - A (1) AYBB {1) = g(Y), 

where Y = X — Xq. From this we see that every solution X of the matrix 
equation ([2|) can be represented in the form (|25|) . 



Based on the following matrix equality 

cf{Y)=g{Y) + {X Q -A^CB^) 
we see that the function g satisfies the condition |T]) iff X Q — A^CB^ 1 ' 

So, the general solution ([23)1 of the matrix equation @ is reproductive iff Xq = 
A^CB^. Therefore, Penrose's general solution (|23|) of the matrix equation 
|2]) is the reproductive solution. 

Remark 2.2. If the condition (16) is not true, the matrix equation (0) is solved 
approximately as described in \2&] . and fH/.O 

2.3. Using the form of the matrix C which is necessary and sufficient for satisfy- 
ing the condition (|17[) we give the necessary and sufficient form for a particular 
solution Xq of the matrix equation @ so that the general solution (l2"5j) of the 
matrix equation @ is reproductive. 
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Theorem 2.4. Let Xq is any particular solution of the matrix equation (0). 
The general solution \25\) of the matrix equation 0) is reproductive iff 



X = Pi 



Ci 


C X Y X 


X2CI 


X2C1Y1 



where P\, Q2, X2, Y\ are the matrices from \10\) and C\ is the submatrix of the 
matrix C and it has the following form: 



cx,x 



C a ,X 



Cx,b 



Ca,b 



where Ci j are some elements o/C. 



Proof. The statement follows from Theorem 



awcbw = Pi 



la 


Xx ' 


x 2 


x 3 



because 
QxCP 2 



' h 


Yx " 




Y 3 . 



(8) 



Pi 



= Pi 

HOI 

= Pi 



la 


Xx ' 


( 






'Pi 


h 


Yx 


x 2 


x 3 




Y 2 


Y 3 


la 


Xx ' 


OxcK 




h 


Yx 




Q: 




x 2 


x 3 


[ ^ 


Y 3 




I 


la 


Xx ' 




r Cx 


( 


" 






h 




Yx 




x 2 


x 3 












[ Y 2 


Y 3 


Q2 



Q2 



Cx 



X2C1 



CxYx 



XiC\Yx 



Q2. ♦ 



Remark 2.3. In the paper fSa l authors proved that there is a particular solution 
Xq of the matrix equation (2) so that Xq ^ A^CB^ for any choice of {1}- 
inverses and .() 



2.4- In this part of the paper we analysed solutions of some matrix systems 
using the concept of reproductivity. 

Application 2.1. In Jj/ R. Penrose studied a matrix system 

(26a) AX = B A (266) XD = E, (26) 

where A, B, D and E are given complex matrices corresponding dimensions. 
He proved that 

Xx=A {1) B + ED {1) -A {1) AED {1) (27) 

is one common solution of the matrix equations (26a) and (26b) provided AE = 
BD and the matrix equations (26a) and (26b) are consistent. 
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In |U/ A. Ben-Israel and T.N.E. Greville proved that the matrix equations (26a) 
and (26b) have a common solution iff each equation separately has a solution 
and AE = BD. Also, they proved that if Xq is any common solution of the 
matrix equations (26a) and (26b), the general solution of the matrix system 
\26)) is given by the formula 

X = g(Y) =X + (I- A^A)Y(I - (28) 

where Y is an arbitrary matrix corresponding dimensions. 

By multiplying the matrix equation (26a) by D on the right and the matrix 
equation (26b) by A on the left we get 

AXD = BD, AXD = AE. (29) 

From this we conclude that AE = BD. 

We will prove that if the matrix system H26)) is consistent, the general reproduc- 
tive solution is given by the formula 

X = f{Y) = A^B + ED& - + (I - A^A)Y(I - DD<V), (30) 

where Y is an arbitrary matrix corresponding dimensions. 

If the matrix system \26]) is consistent, the following equivalence is true. 

{AX = B AXD = E) X = f{X) (31) 

The direct implication of 13 1\) follows by the following implications (see Re- 
mark 2. 1.): 

AX = B => X = f 1 (X) = A w B + X-A {1) AX, 
XD = E =^> X = f 2 (X) = ED {1) +X-XDD ( - 1 ' ) , 
AXD = BD = AE => X = f 3 {X) = A {1) AED {1) +X- A {1) AXDD {1) . 
From the previous implications we get that 

(AX = B A XD = E) => X = f(X) = f 1 (X) + f 2 {X)-h{X). 

The reverse implication of 131)) is trivial. Notice that the function f is repro- 
ductive. Therefore, if the matrix system 1126]) is consistent, it is equivalent to the 
reproductive matrix equation X = f(X). Based on Theorem 1.2. we conclude 
that X = f(Y) is the general reproductive solution of the matrix system \26)) . 
It is easy to see that if there is a particular solution Xq so that Xo ^ X\, then 
X = g(Y) is the general non-reproductive solution. 

Application 2.2. Let A 6 C' iX ™ be a singular matrix. In this section we 
consider a matrix system 

AX A = A A AX = XA. (32) 
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The consistency of the matrix system is determined by Theorem 1 in Urn . 
Let A is commutative {1} -inverse. Based on the reproductivity, we give a new 
proof that the formula 

X = f(Y) = Y + AAA - AAY - YAA + A AY A A, (33) 

where Y is an arbitrary matrix corresponding dimensions, represents the general 
solution of the matrix system h32]) . The formula &33\) is listed by Theorem 3 in 
the paper FZil /. 

// the matrix system \32]) is consistent, the equivalence 

(AX A = A A AX = XA) X = f(X) (34) 

is true. 

The direct implication of \3J$ is based on the following equalities: 

A^A^A = AAXAA = AAX^AA = AXAAA = AXA = AAX 

(=AXA) (= AA) (= X A) (= A) (=AX) 

and 

AaxaA = aAxaA^aaaxA^ ax / A = xaA. 

(= AA) (= AX) (= A) (= XA) 

From this we get that X = X + AAA - AAX + AAXAA - XAA = f(X) . The 



reverse implication of \3J$ is trivial. Notice that the function f is reproduc- 
tive. Therefore, if the matrix system 132\) is consistent, it is equivalent to the 
reproductive matrix equation X — f(X). Based on Theorem 1.2. we conclude 
that X — f(Y) is the general reproductive solution of the matrix system V32\l . If 
Xq is any solution of the matrix system H32\) . the formula 



X = g(Y) = Y + X - AAY - YAA + AAY AA, (35) 



also determines the general solution of the matrix system \32\) because the equal- 
ity g(Y) — f (Y + Xq) is true. If there is a particular solution Xq of the matrix 
system 132]) so that Xq^ AAA, then X = g(Y) is the general non-reproductive 
solution. The previously described process of proving the generality of the solu- 
tions can also be applied to the matrix system which determines k-commutative 
{1} -inverse from flu l. F/ll /. 
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